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Abstract 

We investigate the motion of test particles in the gravitational field of a static naked singu¬ 
larity generated by a mass distribution with quadrupole moment. We use the quadrupole-metric 
(g—metric) which is the simplest generalization of the Schwarzschild metric with a quadrupole 
parameter. We study the influence of the quadrupole on the motion of massive test particles and 
photons and show that the behavior of the geodesics can drastically depend on the values of the 
quadrupole parameter. In particular, we prove explicitly that the perihelion distance depends on 
the value of the quadrupole. Moreover, we show that an accretion disk on the equatorial plane of the 
quadrupole source can be either continuous or discrete, depending on the value of the quadrupole. 
The inner radius of the disk can be used in certain cases to determine the value of the quadrupole 
parameter. The case of a discrete accretion is interpreted as due to the presence of repulsive gravity 
generated by the naked singularity. Radial geodesics are also investigated and compared with the 
Schwarzschild counterparts. 

PACS numbers: 05.70.-a; 02.40.-k 
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I. INTRODUCTION 


The black hole uniqueness theorems state that in Einstein’s general relativity thec^ 
the most general black hole solution in empty space is described by the Kerr metric [l|, 
which represents the exterior gravitational field of a rotating mass m with specific angular 
momentum a = J/m. In Boyer-Lindquist coordinates {t, r, 0,ip), a true curvature singularity 
is determined by the equation r^+a^ cos^ 9 = 0 that corresponds to a ring singularity situated 
on the equatorial plane 6 = 7r/2. The ring singularity is isolated from the exterior spacetime 
by an event horizon situated on a sphere of radius = m + ^/m? — a^. 

In the case that a? > no event horizon exists and the ring singularity becomes 
naked. Different studies show, however, that in realistic situations, where astrophysical 
objects are mostly surrounded by accretion disks, a Kerr naked singularity is an unstable 
configuration that rapidly decays into a Kerr black hole. Moreover, it now seems established 
that in generic situations a gravitational collapse cannot lead to the formation of a final 
configuration corresponding to a Kerr naked singularity. These results seem to indicate that 
rotating Kerr naked singularities cannot be very common objects in nature. 

The above results seem to corroborate the validity of the cosmic censorship hypothesis 
according to which a physically realistic gravitational collapse, which evolves from a regular 
initial state, can never lead to the formation of a naked singularity; that is, all singularities 
formed as the result of a realistic collapse should always be enclosed within an event horizon; 
hence, the singularities are always invisible to observers situated outside the horizon. Many 
attempts have been made to prove this conjecture with the same mathematical rigor used 
to show the inevitability of singularities in general relativity jo]. No general proof has 
been formulated so far. Instead, particular scenarios of gravitational collapses have been 
investigated some of which indeed corroborate the correctness of the conjecture. 


Nevertheless, other studies indicate that under certain circumstances naked singular¬ 
ities can appear as the result of a realistic gravitational collapse. Indeed, it turns out that 
in the case of the collapse of an inhomogeneous matter distribution, there exists a critical 
degree of inhomogeneity below which black holes form. Naked singularities appear if the 
degree of inhomogeneity is higher than the critical value. The speed of the collapse and the 
shape of the collapsing object are also factors that play an important role in the determina¬ 
tion of the final configuration. It turns out that naked singularities form more frequently if 
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the collapse occurs very rapidly and if the object is not spherically symmetric. 

In view of this situation it seems reasonable to investigate the effects of naked singularities 
on the surrounding spacetime. This is the main aim of the present work. We will study 
a naked singularity without black hole counterpart. In fact, we investigate the quadrupole 
metric (g—metric) which is the simplest generalization of the Schwarzschild metric containing 
a naked singularity. We will see that, starting from the Schwarzschild metric, the Zipoy- 
Voorhees |8|, l9| transformation can be used to generate a static axisymmetric spacetime 
which describes the held of a mass with a particular quadrupole moment. For any values 
of the quadrupole, the spacetime is characterized by the presence of naked singularities 
situated at a hnite distance from the origin of coordinates. In this work, we are interested 
in analyzing the spacetime outside the outer naked singularity. 

We perform an analysis of the motion of test particles (massive and non-massive) outside 
the naked singularity, comparing in all the cases our results with the corresponding situation 
in the Schwarzschild spacetime in order to establish the exact inhuence of the quadrupole 
on the parameters of the trajectories. 

This paper is organized as follows. In Sec. m we present the g—metric and investigate 
its main physical properties. It is shown that it can be used to describe the exterior grav¬ 
itational held of a mass distribution with quadrupole moment. Moreover, we prove that 
the corresponding spacetime is characterized by the presence of naked singularities. In Sec. 
HU we derive the geodesic equations in general and study numerically the conditions un¬ 
der which motion is allowed. Sec. lYl is devoted to the study of circular and non-circular 
equatorial geodesics. We show that the geometric conhguration of an accretion disk located 
around a naked singularity depends explicitly on the value of the quadrupole in a such a 
way that it is always possible to distinguish between a Schwarzschild black hole and a naked 
singularity. This effect was found previously in the held of naked singularities with electric 
charge 10, lll[. Furthermore, in Sec. IVj we study the inhuence of the quadrupole on the 
behavior of radial timelike and null geodesics. Finally, in Sec. IVTl we discuss our results 
and comment on possible future works. 
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II. THE g-METRIC 


In 1917 [1^, it was shown that the most general static axisymmetric asymptotically flat 
solntion of Einstein’s vacnum equations is represented by the Weyl class. In terms of multi¬ 
pole moments, the simplest static solution contained in the Weyl class is the Schwarzschild 
metric which is the only one that possesses a mass monopole moment only. From a physical 
point of view, the next interesting solution should describe the exterior held of a mass with 
quadrupole moment. In this case, it is possible to And a large number of exact solutions 
with the same quadrupole (see 13| and the references cited therein) that differ only in the 
set of higher multipoles. A common characteristic of the solutions with quadrupole is that 


their explicit form is rather cumbersome, making them difficult to 
[l^ . An alternative exact solution was presented by one of us in 


De handled analytically 


ISj by applying on the 


Schwarzschild metric a Zipoy-Voorhees transformation with parameter 6 = 1 + q, where q 
represents the quadrupole parameter. In spherical coordinates, the resulting metric can be 
written in a compact and simple form as 


ds^ = 


2m \ 

1 - —) dt^ 
r J 


r J 


1 + 


m2 sin2 0\V dY 


— 2 mr 


_ 2m 


-I- YdO'^ + Y sin^ OdY 


•( 1 ) 


In the literature, this solution is known as the metric or as the 7 —metric for notational 
reasons [^. We propose to use the term quadrupole metric (g—metric) to emphasize the 
role of the parameter q which determines the quadrupole moment, as we will see below. 

The g—metric is an axially symmetric exact vacuum solution that reduces to the spher¬ 
ically symmetric Schwarzschild metric only for g —)• 0. It is asymptotically flat for any 
finite values of the parameters m and g. Moreover, in the limiting case m —)■ 0 it can be 
shown that, independently of the value of g, there exists a coordinate transformation that 
transforms the resulting metric into the Minkowski solution. This last property is important 
from a physical point of view because it means that the parameter g is related to a genuine 
mass distribution. 

An important quantity that characterizes the physical properties of any exact solution is 
the Arnowitt-Deser-Misner (ADM) mass. If one tries to use the common formula [^, ^ 
for calculating the ADM mass, one would obtain incorrect results because these formula are 
adapted to a particular coordinate system. Therefore, we perform here a detailed calculation 
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of the ADM mass, using the original approach. All the details are given in the Appendix. 
In the case of the g—metric o. the hnal result is simply 

Madm = m{l + q) . (2) 


It follows that if we take the parameter m as positive, the parameter q must satisfy the 
condition g > — 1 in order for the ADM mass to be positive. Nevertheless, one can choose 
a negative m so that for g < — 1 the mass is still positive. For the sake of simplicity we 
assume from now on that m is positive. 

We also calculate the multipole moments of the g—metric by using the invariant dehnition 
proposed by Geroch jl^. The lowest mass multipole moments M„, n = 0,1,... are given by 


m 


Mo = (l + g)m, Ma = —^g(l + g )(2 + g) , 


( 3 ) 


whereas higher moments are proportional to mg and can be completely rewritten in terms 
of Mq and M 2 . This means that the arbitrary parameters m and g determine the mass 
and quadrupole which are the only independent multipole moments of the solution. In the 
limiting case g = 0 only the monopole Mq = m survives, as in the Schwarzschild spacetime. 
In the limit m —)■ 0, with g 7 ^ 0, and g —>■ —1, with m 7 ^ 0, all multipoles vanish identically, 
implying that no mass distribution is present and the spacetime must be flat. Furthermore, 
notice that all odd multipole moments are zero because the solution possesses an additional 
reflection symmetry with respect to the equatorial plane 6 = tt/2. 

We conclude that the above metric describes the exterior gravitational held of a static 
deformed mass. The deformation is described by the quadrupole moment M 2 which is 
positive for a prolate mass distribution and negative for an oblate one. Notice that the 
condition g > — 1 must be satished in order to avoid the appearance of a negative total mass 
Mq. Therefore, in the interval g G (—1, 0) the g—metric describes a prolate mass distribution 
and in the interval ( 0 , cxo) an oblate one. 

To investigate the structure of possible curvature singularities, we consider the 
Kretschmann scalar K = . A straightforward computation leads to 

16m^(l + g)^ (r^ — 2mr + rn? sin^ 6 i) 2 ( 2 q+A)-i 


K = 


^4(2+2g+g2) 


(1 — 2m/r) 2 ( 9 ^+9+1) 


-L(r, 0 ) , 


( 4 ) 


with 


L(r, 6) = 3(r — 2m — qmY{r'^ — 2mr + m^ sin^ 9) 

+g(2 + q)m? sin^ 9[q{2 + g)m^ + 3(r — m)(r — 2m — qm)] 


( 5 ) 
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FIG. 1: Curvature singularities of the g—metric. The outer singularity at r = 2m and the singu¬ 
larity at the origin r = 0 exist for any value of the parameter q. The singularity which depends on 
the angle 6 exists only for certain values of the quadrupole (see text). 

In the limiting case g = 0, we obtain the Schwarzschild value K = /r^ with the only 

singularity situated at the origin of coordinates r —)■ 0. In general, we can see that the 
singularity at the origin, which occurs at 7-4(2+2i}+(j2) _ jg pj-gggni; for all real values of q. 
Moreover, an additional singularity appears at the radius r = 2m which, according to the 
form of the metric o. is also a horizon in the sense that the norm of the timelike Killing 
vector vanishes at that radius. Outside the hypersurface r = 2m no additional horizon 
exists, indicating that the singularities situated at the origin and at r = 2m are naked. In 
addition, a curvature singularity occurs at the surface determined by the equation 

— 2mr -|- sin^ 0 = 0 (6) 

under the condition that the value of the quadrupole parameter is within the interval 
q G (-1,-1 -f ■\/3/2]\{0}. We conclude that the main consequence of the existence of 
a quadrupole determined by the parameter q is that the Schwarzschild horizon becomes a 
naked singularity. The geometric structure of the curvature singularities of the g—metric is 
illustrated in Fig. [H 

Several physical properties of the g—metric have been investigated in the literature (see. 
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for instance 


, Q, Q, £- 


26| and references therein). In this work, we are interested in 


continning the analysis of the corresponding gravitational held as observed by exterior test 
particles. 


III. GEODESIC MOTION 


Consider the trajectory x“(r) of a test particle with 4-velocity /dr = x". The 

moment p" = /rx" of the particle can be normalized so that 


= e, 


( 7 ) 


where e = 0,1, —1 for nnll, timelike, and spacelike cnrves, respectively 10|. For the g—metric 
we obtain from ([7]) that 


1 + 


2 ■ 2 n \ - 9 ( 2 + 9 ) 

m sm 6 ^ ' 


— 2mr 


= E^ - 4)2, 


where 


2m \ 

4)^ = ( 1 - 

r / 


1+9 


2, 2 m\-V 

^ ( 7 -) ( 7 “I—2 o ) ^ —2 ■ 2 a 

r J \ — 2mr J sm 6 


( 8 ) 


2 m \ 

1 -I + e 

r / 


( 9 ) 


and we have used the expression for the energy E = fiE and the angular moment I = fil of 
the test particle which are constants of motion 

1+9 


E = ^ ~ ?") 


( 10 ) 

( 11 ) 


associated with the Killing vector helds = dt and respectively. For the sake of 

simplicity we set p = 1 so that E = E and 1 = 1. 

In addition, the equation for the acceleration along the polar angle can be expressed as 

9(2+9) 


0 _ _ ( q{2+q)'m? sing cos 6> \ f _i 

Y r‘^—2m,r+m? sin*^ d j yr'^—2mr J ' 


—2mr+m^ sin^ 9 
r^— 2 mr 


sin 9 cos 9 


•2 


r® —(4+g)mr^ + (2(2+9)+(l+g)^ sin^ (l+9)(2+9)m® sin^ 0 . ^ 

(r^—2mr)(r^—2mr+m^ sin^ 0 ) 


( 12 ) 


From this equation it follows that the acceleration 9 depends on the polar angle. Consider 
for instance the motion around the central object with initial values r = 9 = 0 and 7 ^ 0. 
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Then, Eq. dT^ reduces to 


/ - 2 Q \ -9(2+9) 

( 1 +-r-] 0 — sin 6* cos 00^ = 0 , (13) 

y 7-2 _ 2mr J 

It then follows that the only places at which the acceleration vanishes are 0 = 0 and 9 = 7r/2, 
showing that the polar and equatorial planes are geodesic planes at which the investigation 
of the test particle motion can be considerably simplihed due to the symmetry of the conhg- 
uration. For all the remaining values of 0, there is always a non-zero value of the acceleration 
that induces a velocity in the direction of 0 towards the equatorial plane or the poles. 

Equation (jH]) can be considered as describing the motion along the radial coordinate in 
terms of the “effective potential” However, the velocity 0 enters the expression for 
explicitly so that strictly speaking it is not a potential. Nevertheless, its dependence on the 
coordinates can shed some light on the behavior of the geodesics for particular values of 0. 
In Fig. [21 we plot the behavior of for a particular initial velocity 0 = 1 and different 
values of the q parameter. We see that at the singularity r = 2m the behavior of the 
potential drastically depends on the value of q. Indeed, for positive values of q the potential 
either tends to inhnity near the singularity or crosses it and becomes divergent inside the 
singularity. Instead, if q is negative, the potential tends to inhnity as the singularity is 
approached from outside. In the case of vanishing quadrupole, the singularity turns into 
a horizon and the potential crosses it until it becomes inhnity near the central singularity 
located at the origin of coordinates. 

An interesting particular case of the geodesic motion is that of circular orbits and their 
stability. The geodesic equation in this case is given as in Eq.([8]) with r = 0, i.e. 
whose solutions strongly depend on the value of the velocity 0. In fact, the condition for 
circular orbits can be expressed as 

P / 2my+^'^ 2my+‘'' 

sin^ 0 \ r J \ r J 

(14) 

We perform a numerical investigation of this condition for a particular radius and diherent 
values of the quadrupole parameter in Fig. O If we consider, for instance, the case of a 
circular orbit with radius r = 2.35m, we can see that such a motion is not possible for 
arbitrary values of q and 0. In fact, for most negative values of q no circular orbit with this 
radius is allowed, because 0^ is negative for all values of 0. Instead, for all positive values 


0 ^ = {m—2mr) I 1 -|- 


,2 Q \ 9(2+9) 


m sm' 


2mr 
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FIG. 2: Behavior of the “effective potential” in terms of the radial distance r for the particular 
velocity 6 = 1 with 0 = 7r/4, I = 15 and different values of q. The effective potential for the 
Schwarzschild solution g = 0 on a fixed plane with 0 = 0 is also depicted for comparison. 
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FIG. 3: Angular velocity 9 in terms of the angle 6 G (0, vr) for different quadrupoles. Here, we 
consider timelike geodesics (e = 1) with m = 1, I = 5. For concreteness, the radius of the orbit has 
been set to r = 2.35. 
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of q it is always possible to find an interval of 6 within which it is possible to have circular 
orbits. For other values of the orbit radius similar results can be found. 


The study of the stability of circular orbits is important to establish the possibility of 
having accretion disks around the central gravitational source. In particular, the radius of the 
last stable circular orbit determines the inner radius of the accretion disk, and corresponds 
to an inflection point of the potential, i.e., the point where the conditions /dr = 0 and 
52 <J) 2/^^2 _ g satished. We use these two conditions to hnd the explicit values of P and 
6*^ for the last stable circular orbit and obtain 


em(l + q) sin^ 6 
r — (3 + 2q)m 


2m\ ^ 

^ {r — m) G 

V r ) 

r(r — 2m) [r — (3 + 2g)m] H 


(15) 


and 

.2 em(l + g)(l — — 2mr + m? sin^ 

— (3 + 2g)m](r^ — — 2mr + m‘^{l + qY sin^ 0] H' 

where 

G = r^ — 2m{4 + 3q)r + 2rrY{2 + g)(3 + 2q) , 


(16) 


(17) 


LT — -I j_ r—m _ 

' (r^ —2mr) [r— (3+2g')m)] 



6m(3 + 2q)r + 4m^(2 + g)(3 + 2q) 


I 2g(2+g)(l+g)^m'^ sin^ 6{r—m)[r—{3-\-2q)m] 
(r^—2mr+m^ sin^ 0)[r^—2mr+m^(l+g)^ sin^ 0] 


(18) 


In the limiting case of the Schwarzschild spacetime (g = 0), due to the spherical symmetry 
we can set 9isco = 0 and so we obtain the value rf/Y, = 6m, as expected. The value of 
determines the angular momentum of the test particle on the last stable circular orbit 
and the velocity of the orbit with respect to the polar angle. In Fig. 01 we analyze 
numerically the above equations for different values of the quadrupole. The condition for the 
existence of a last stable circular orbits is that both and 9f^^^ be positive and hnite. From 
Fig. 0]we see that for q = 0.5 and q = —0.5 these conditions are not satished. For q = 0.1 
and q = —0.1, however, there exists an interval of r in which both and 9f^^^ are positive. 
The hrst thing one can notice is that this interval is always inside the Schwarzschild radius 
rffY, = 6m. This means that the quadrupole moment diminishes the value of the radius of 
the last stable circular orbit. For a particular positive value of 9f^^^ one can hnd from the 
graph the corresponding radius risco- If this value happens to correspond to a positive value 
of then there exists a stable circular orbit with that radius and that angular momentum. 
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FIG. 4: Angular momentum and polar angle velocity of the last stable circular timelike orbit as a 
function of the radial coordinate for different values of the quadrupole parameter q. Here we set 
9 = IT/A. The value of has been rescaled in all the graphs for comparison with the graph of 


The motion of test particles for arbitrary values of 6 and 6 can be investigated by per¬ 
forming a numerical integration of the geodesic equations l|5 |) - ffT^ . Several cases must be 
considered, depending on the initial velocity 6 and the initial plane 6. This requires a de¬ 
tailed numerical analysis in which several aspects must be considered like the stability of the 
trajectories and the appearance of chaotic motion. We will present these results elsewhere. 
In this work, we will focus on the analytic investigation of the most important families of 
geodesics of the q— metric. 
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IV. EQUATORIAL GEODESICS 


As mentioned above the equatorial plane 0 = 7r/2 is a geodesic plane due to the symmetry 
of the gravitational source. In this case, the geodesic equations reduce to 


t = E[l - 

r / 


(19) 




r J 


1 + 


m 


2mr 


-q{2+q) 


r2 = I l_ 


2m 


9+1 r 




^ 2m V 

1- 1 +e 

r 


( 20 ) 


( 21 ) 


To investigate the behavior of the geodesics it is convenient to introduce the new coordinate 
u = 1/r and use the azimuthal angle ip instead of the affine parameter r so that 


^ = lu^l-2mur^ . 
dr dp 

Then, the geodesic equation fl2lD can be expressed as 

(1 - 7nufi<3+i) [iP) ^ *“’'**’ 


( 22 ) 


(23) 


where 

F{u, <?) = ^ “ (1 “ 2muY~^'^ — 2muY + e] } . (24) 

The properties of the geodesic motion are then determined by the function F{u,q) which 
must be positive for the velocity to be a real valued function. Furthermore, the zeros of 
F{u,q) are the points where the radial velocity vanishes, i.e., the perihelion or aphelion 
distance in the case of bounded orbits. In turn, the behavior of F{u, q) can be explored by 
comparison with the function f{u) which is dehned as 


/(«) = F(a,q = 0) = p - (1 - 2rrm)(l^ii^ + ()] 


(25) 


This means that the function f(u) determines the behavior of the geodesics in the case 
of the Schwarzschild spacetime and, consequently, a comparison between F{u,q) and f{u) 
will determine the main differences between geodesics of the g—metric and those of the 
Schwarzschild spacetime. It is easy to show that 


F{u,q)>f{u) for g>0 and F{u,q) < f{u) for g<0 (26) 
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FIG. 5: Behavior of the function F{u,q) for particular values of the parameters m = 1, E = 1, 
and I = 15, and e = 1. 


in the interval u G (O, i.e., in the spatial region located between the outer singularity 

(r = 2m) and infinity. This behavior is illustrated in Fig. [5l This simple observation leads 
to an interesting physical result. Suppose that the zeros of f{u) are ui and U 2 with ui < U 2 
so that the trajectory is bounded between ri = 1/ui and r 2 = l/rt 2 ; consequently, the 
perihelion distance is r 2 < ri. Suppose that the value of q is so chosen that F{u, q) has also 
two zeros at uj < u' 2 - Then, from the behavior of F{u, q) (cf. Fig. [5]) we conclude that 

Ml < u'l , i. e. ri > rj and < U 2 , i. e. > r 2 for g > 0 . (27) 


This means that the perihelion distance for positive q is greater that the one for vanishing 
quadrupole. A similar analysis shows that the perihelion distance decreases for negative 
values of the quadrupole. We conclude that the presence of the quadrupole in the g—metric 
leads to a change of the perihelion distance, an effect that has been predicted also for other 
metrics with quadrupole moment [l4, 27]. Some examples of geodesics for different values of 
the quadrupole parameter are given in Figs. [MHl where for simplicity the radial coordinate 
has been redehned as r/m for all the trajectories. 

In Fig. [6], we consider unbounded Schwarzschild orbits with non-vanishing initial radial 
velocities under the influence of the quadrupole. We see that for the chosen initial radial 
and angular velocities all the test particles are able to escape from the gravitational held 
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FIG. 6: Influence of the quadrupole on unbounded Schwarzschild orbits with non-vanishing initial 
radial velocity (r(0) 7^ 0). The initial point <^(0) = 0 and r(0) = 20 is the same for both plots. 
Initial velocities: f(0) = —2 and (/9(0) = 0.04305 (left plot); r(0) = —2.5 and :/9(0) = 0.039 (right 
plot). 



FIG. 7: Influence of the quadrupole on unbounded Schwarzschild orbits with vanishing initial radial 
velocity (r(0) = 0). Initial conditions: (/p( 0) = 7r/2, r(0) = 7, (^(0) = 0.07145 (left plot); (/?(0) = 0, 
r(0) = 9, <^(0) = 0.075 (center plot); (/^(O) = 7r/2, r(0) = 9, 99(0) = 0.07145 (right plot). 

of the naked singularity. The value of the quadrupole determines only the direction along 
which the particle escapes towards inhnity. 

In Fig. [71 we investigate how the quadrupole acts on unbounded Schwarzschild orbits 
with zero initial radial velocity. First, we study the change of a stable circular Schwarzschild 


14 























30 



FIG. 8: Influence of the quadrupole on Schwarzschild bounded orbits with vanishing initial radial 
velocity (f(0) = 0). The initial conditions are (/?(0) = 0, r(0) = 7, and <^(0) = 0.08 for all the 
trajectories. 

geodesic under the influence of a quadrupole (left panel). For a small quadrupole (g = 0.009), 
the geodesic moves slowly towards the central singularity and Anally reaches a stable circular 
orbit with a radius which is smaller than the initial Schwarzschild radius. In the central 
and right plots, we consider the same values for the quadrupole moment, but diAerent 
initial angular velocities 0(0). In both cases, the test particles move along bounded orbits 
with diAerent perihelion shifts. From the behavior of these trajectories, we conclude that 
the quadrupole is able to transform unbounded Schwarzschild trajectories into bounded 
trajectories. 

In Fig. IH], we consider a Schwarzschild bounded orbit under the inAuence of the 
quadrupole with zero initial radial velocity and the same non-zero value for the initial angu¬ 
lar velocity. The left panel shows the Schwarzschild geodesic, whereas the central and right 
plots are for a particular negative and positive value of the quadrupole, respectively. In 
this case, we see that the quadrupole does not aAect the bounded character of the geodesic. 
However, it can drastically modify the geometric structure of the trajectory. 

In general, we see that the quadrupole of the naked singularity always aAects the motion 
of test particles. The speciAc change can manifest itself in diAerent ways, depending on the 
explicit value of the quadrupole. 
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FIG. 9: The effective potential for timelike circular geodesics on the equatorial plane as a function 
of the radius for different values of the quadrupole parameter. Here we set = 30 for concreteness. 

A. Circular orbits 


Circular orbits on the equatorial plane can be investigated analytically. In fact, in this 
case the motion under the condition r = 0 is equivalent to the motion of a test particle in 
the effective potential [cf. Eq. fl2lD ] 



2m 


r / 


p 



(28) 


The presence of the quadrupole parameter influences the behavior of the effective potential, 
as can be seen in Fig. |9l The effective potential of the Schwarzschild is also shown for 
comparison. For positive values of q the value of the effective potential at a given point 
outside the outer singularity is always less than the Schwarzschild value. For negative values 
of the quadrupole the opposite is true. This indicates that the distribution of circular orbits 
on the equator of the g—metric can be modihed drastically by the quadrupole. We will now 
explore this point in detail. 

The explicit value of the angular momentum of a test particle along a circular orbit can 
be derived from the condition dV^j^/dr = 0. A straightforward computation yields that this 
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condition is satisfied if the angular momentum is given by 


„ m(g + l)(l-gf) 
r — (3 + 2q)m 


( 29 ) 


an expression from which it follows immediately that the radius of any circular orbit must 
satisfy the condition 

Tc > m{3 + 2q) . (30) 


This radius is greater than the Schwarzschild radius for a positive quadrupole parameter. 
However, for negative values of q the radius is smaller than the Schwarzschild value and, in 
principle, can be made as close as possible to the outer singularity located at r = 2m. Below 
we will show that in the limiting case Tc = m(3 + 2q), the orbit becomes lightlike. 
Furthermore, the energy of a test particle in circular motion can be expressed as 




2m m(g + 1) 

r J r — m(3 + 2q) 


(31) 


which is positive only within the range of allowed radii fl30|l . To completely characterize the 
parameters of circular orbits, we also calculate their angular velocity Q{r) and period T(r) 


and obtain 


T(r) 


respectively. 


f2(r) = ip = - 
r 


m(l + 9) (1 - ^-f)" 

r — m(3 + 2q) 


1/2 


—dip 

0 



27rr^/2 


r — m(2 + q) 
m(l + q){r — 2m) 




(32) 

(33) 


B. Stability analysis 

We now analyze the stability properties of the circular motion. In particular the last 
stable circular orbit is determined by the inflection points of the effective potential (|28|) . i.e., 
by the zeros of the equation 

m{l + q) (r^ — 8 mr — 6 mqr + 4 m^g^ + 14 m^g + 12 m^) 
dr'i V / r‘^\r — m{3 + 2g)] 

where we have replaced the value for the angular momentum of the circular orbit fl2^ . Then, 
we obtain the radii 

hIo = "^(4 + 3g ± v/5g2 + lOg + 4) (35) 
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FIG. 10: The last stable circular radius for a timelike particle as a function of the quadrupole 
parameter q. The critical radius = m(3 + 2q) and the outer singularity = 2m are also 
plotted. The left plot is a zoom of the region of intersection. 

which are positive and lie outside the outer singularity for all allowed values of q. The 
behavior of is plotted in Fig. [10] where we also plot the radius rc that denotes the 
minimum value at which the angular momentum is positive. It follows from the graph 
that it is necessary to consider three different intervals: I for q G (oo,—0.5], II for q G 
(—0.5, —1 + l/\/5] and III for q G (—1 + l/\/5, —1). The particular values of q used to 
determine the different intervals follow from the conditions r_ = {q = —0.5) and r+ = r_ 
(q = -l + l/x/5^ -0.5527). 

In the interval I, the only valid radius is which is greater (smaller) than the 

Schwarzschild radius for positive (negative) values of the quadrupole parameter q. No¬ 
tice that on the boundary of this interval (g = —0.5) the zero at cannot be considered 
because it is located on the outer singularity. In this interval we have that > 0 

so that all the circular orbits are stable. If we imagine a hypothetical accretion disk made 
of test particles around the central source so that the inner radius of the disk coincides with 
^Zsco! quadrupole in the interval I consists in changing the inner radius. From 

an observational point of view, this implies that by measuring the inner radius of the disk, 
one can determine the value of the quadrupole. The smallest disk inner radius corresponds 
to the value ~ —0.5) = 3m whereas the outer radius can in principle be extended to 

infinity. 

In the interval II with q G (—0.5, —1 -f l/\/5), the second derivative of 17^^ is negative 
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FIG. 11; The angular momentum of test particles around a naked singularity with q = —0.52. The 
left plot is a zoom of the region where the internal stable ring is located. 

within the two zeros located at and This means that circular orbits are unstable 
within this interval. The radial extension L of the unstable region is 

, , - m for q = —0.5 

T = Gtco - G.co = 2m + lOg + 4 = <^ (36) 

[O forg = -l + l/v^ 

The interesting fact in this case is that there exists an additional stable region located be¬ 
tween the outer singularity and (cf. Fig. [10]) which is separated by the unstable region 
(^Zsco) ^itco) from the exterior stable region with r > An accretion disk around such 

an object would consist of an external disk which can extend from to infinity and an 
internal ring located inside the region (2m, m-|-3m/\/5 ~ 2.34m), where the outer boundary 
corresponds to the point where Since the internal stable ring is located so close 

to the outer singularity one could imagine that the angular momentum of the test particles 
should be very high in order to maintain the orbit. A detailed analysis, however, shows that 
this is not true. In Fig. [TTl we plot the value of the angular momentum as a function of 
the radial distance for a particular quadrupole parameter which allows the existence of an 
internal ring. We can see that the maximum value of the angular momentum is reached on 
the external boundary of the ring and then it decreases as the singularity is approached. 
The explicit value of the angular momentum is not high and in fact it is comparable with 
the values for test particles located inside the external stable disk at distances of several 
times m. The peculiarity of the angular momentum is that its value decreases as the sin¬ 
gularity is approached, an effect that contradicts the physical expectations for an attractive 
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gravitational field. The only possible explanation for this unusual behavior is that there is 
an additional force that compensates the gravitational attraction. Similar effects has been 
found in the gravitational held of other naked singularities and have been interpreted as a 


manifestation of repulsive gravity 


3,|nl,£, 


29| . An alternative approach in w 


30 


rich repulsive 


3ll. 


gravity is dehned in terms of curvature invariants leads to similar results 

In the interval III, with q G (—1 + l/\/5, —1), there are no last stable orbits. All the 
circular trajectories are allowed, in principle, starting at the outer singularity r = 2m. 
An analysis of the angular momentum shows that it diminishes as the naked singularity is 
approached. Again, the simplest explanation is to assume the presence of repulsive gravity. 
An accretion disk around this type of naked singularities would have a continuous structure 
which extends from the singularity to inhnity. 


C. Circular null geodesics 


In the case of null geodesics (e = 0), the effective potential (1281) reduces to 






2m 
r / 


(37) 


and its behavior for a given P is similar to that shown in Fig. [91 The important quantities 
for the analysis of circular motion of photons are 


91/2 2 / 2 / 2mY\ 

= [r-m{3 + 2q)], 


(38) 


and 


9r2 



2m 
r / 


[3r^ — 12qmr — 18mr + 2%qmP + 8q‘^m^ + 24m^] . 


(39) 


The hrst derivative vanishes for r = r.y = m(3 + 2q), an expression which reduces to the 
standard Schwarzschild value for g = 0, as expected. We see that the circular orbit located 
at r.y = m(3 + 2q) corresponds to the trajectory of a photon. This is also the limiting radius 
for timelike circular geodesics, as explained above. 

By inserting this radius value into the second derivative of the effective potential, we 
obtain an expression which is negative for all values of q. This implies that the circular 
orbit with radius = m(3 + 2q) is unstable, independently of the value of the angular 
momentum. In contrast to the Schwarzschild spacetime, where there is only one circular 
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orbit with = 3m, in the case of the g—metric the radius = m(3 + 2q) can take any 
value within the interval {2m, oo), depending on the value of the quadrupole parameter. 
Nevertheless, for fixed values of m and q only one circular orbit is allowed. 


V. RADIAL GEODESICS 


In the Schwarzschild spacetime, radial geodesics representing the free fall of test particles 
are straight lines that connect the initial point in spacetime with the singularity located at 
the origin of coordinates. Therefore, it is useful to study radial geodesics in the spacetime 
described by the g—metric in order to evaluate the influence of the quadrupole. We will 
consider the free fall of test particles with vanishing angular momentum (/ = 0). Then, the 
motion along the radial coordinate is governed by the equation 

m^ sin^ 9 


r = - 1 + 


|(2+g) 


- $2 , 


— 2mr ^ 

where the minus sign indicates that the particle falls inward and 


(40) 


r J 


^ 2m\ 

r J 


-Q 


1 + 


2 ■ 2 n \ -'j(2+g) 

m sm 9 \ 


2mr 


r+ e 


(41) 


In addition, we must take into account that an arbitrary polar plane {9 =const) is not 
necessarily a geodesic plane. This means that the radial geodesic is determined by Eq. fHOj) 
together with the equation flT^ for the evolution along the angle 9. We illustrate the 
result of the integration of these two differential equations in Fig. [121 The effect of the 
quadrupole parameter g becomes plausible from the behavior of the geodesics. For sources 
with negative quadrupole parameter (prolate sources), the geodesics deviate from the initial 
angle 9 = ^ towards the axis of symmetry. In the case of positive quadrupole g (oblate 
body), the deviation is towards the equatorial plane of source. This result is in accordance 
with our physical expectations. Indeed, an oblate source is larger on the equatorial plane 
and so one would expect that the mass density is larger on the equator, generating a larger 
gravitational interaction. In the case of prolate bodies, the same intuitive reasoning is 
valid for the symmetry axis. This result reinforces our interpretation of the parameter g as 
determining the quadrupole moment of the gravitational source. 

In the case of radial geodesics, it is interesting to compute the coordinate time, t, and 
the proper time, r, that passes along the trajectory of a test mass between two radii rq and 
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FIG. 12: Free fall of test particles with initial point located at r = 4 (m = 1) and 6 = ^. The 
singularity at r = 2m and the radial geodesic for the Schwarzschild spacetime are also shown for 
comparison. 


Tf. For simplicity, we consider the case of radial geodesics on the eqnatorial plane where the 
geodesic eqnations rednce to 


r" = 1 + 


m 


2 \ <?(2+g) 


— 2m 


) '■-( 


r ) 


1+9 


t = E 1- 


2m 


-l-q 


(42) 


It is then easy to obtain the relationships 


At = 


fTf 


r ) 


9+1 


1 - 1/2 


1 + 


m 


— 2mr 


-59(2+9) 


dr , 


(43) 
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and 


At = 


In 


r / 


r / 


- 1/2 


1 + 


m 


— 2mr 


~5'?(2+9) 


dr , (44) 


for the proper time and coordinate time, respectively. Let ns consider the case of a particle 
located initially at r* —)■ cxd so that = 1. If the hnal state is at the singnlarity Vf = 2m, 
we obtain 

/'2m\ "5^+9) 


^2m — 


[-) 


1 + 


m 


dr . 


(45) 


— 2mr ^ 

It is not possible to hnd an analytical expression for this integral. Nevertheless, one can 
evalnate it either nnmerically or by means of a Taylor expansion aronnd the valne q = 0 and 
r = 2m. As a resnlt we obtain that the proper time is hnite for all the allowed valnes of q. 
The same resnlt is obtained for the Schwarzschild spacetime (g = 0). This results indicates 
that the quadrupole does not affect the finiteness of the proper time that is necessary to 
reach the hypersurface r = 2m. 

In a similar manner, for the coordinate time we obtain 

- 1/2 


(*2m 


tom. — E 


E^ -\l- 


2m\ 

r J 


1+9 


^ 2m\ 

r J 


-i+V 


2m m" 
1 -+ ^ 


2 \ - 1 ( 2 + 9 ) 


(46) 


Again, a Taylor expansion of this quantity reveals that it is hnite for g < 0 and inhnite for 
g > 0. This means that a positive quadrupole does not change the main property of the 
coordinate time of the Schwarzschild spacetime. However, in the case of negative quadrupole 
the coordinate time completely changes so that an observer at inhnity would observe how 
the test particle reaches the outer singularity located at r = 2m. 

Let us now consider the radial motion of photons. Suppose that a photon is emitted at a 
radius rem with frequency z/gm and received at a radius r^ec with frequency Urec- Then, the 
redshit z in a static spacetime is determined by the relationship [l^ 


z + 1 = 




9tt\ 


9 tt\e 


1/2 


1 — 2m/rr 
1 — 2m/re 


5(1+9) 


(47) 


If the emission occurs near the outer singularity {rem ~ 2m), the redshift diverges, indepen¬ 
dently of the value of g. This means that for an observer located outside the radius r = 2m, 
it is not possible to receive information from the singularity. In this sense, the singularity 
remains invisible for external observers. 
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VI. CONCLUSIONS 


In this work, we investigated the motion of test particles in the gravitational held de¬ 
scribed by the g—metric, the simplest generalization of the Schwarzschild metric which 
contains a quadrupole parameter. It was shown that the g—spacetime possesses two singu¬ 
larities at r = 0 and r = 2m which are not covered by a horizon, i.e., they are not isolated 
from the exterior spacetime. For certain values of the parameter g, a third dependent 
singularity appears inside the spatial region contained within the radial interval (0,2m). If 
we limit ourselves to the region with r > 2m, the geodesic equations determine the motion of 
test particles around the naked singularity located at r = 2m. First, we studied the bounded 
motion on an arbitrary polar plane 6 ^ ^■,t^/2, and showed that in general the plane of the 
orbit moves towards the equator or to the axis of symmetry, depending on the value of the 
quadrupole. We performed a brief analysis of the conditions under which circular orbits are 
allowed for arbitrary values of the polar angle 9, the angular momentum of the particle and 
the quadrupole parameter. 

We studied in detail the motion on the equatorial plane 9 = 7r/2. We integrated numeri¬ 
cally the geodesics equations and found bounded and unbounded orbits for different values 
of the quadrupole. It was also shown explicitly that the presence of the quadrupole leads 
to a change of the perihelion distance of a bounded orbit. This effect could in principle be 
used to determine the quadrupole moment of the gravitational source. 

The analysis of circular orbits and their stability leads to a classihcation of naked singular¬ 
ities according to which in certain cases it is possible to determine the quadrupole parameter 
by measuring the inner radius of the accretion disks made of test particles only. Suppose, 
for instance, that an accretion disk is detected with inner radius Tinner within the interval 
(3m, cx)). Then, the naked singularity is of type III with a quadrupole parameter that can 
be determined numerically by using the formula Tinner = m(4 -|- 3g -|- A/Sg^ -|- lOg -|- 4) and 
the result must be a value contained within the interval g G (—0.5, cxd). If the value of the 
accretion disk inner radius is inside the interval (2.34m, 3m), the naked singularity belongs 
to the class II, and the quadrupole parameter can be determined again from the formula 
Tinner = TTi^A -|- 3g -f + lOg -f 4). lu this case, the quadrupole parameter must be con¬ 
tained within the interval (-0.5527,-0.5), and there must exist an additional inner ring 
located inside the spatial region (2m, 2.34m), whose exterior radius can be determined by 
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the formula r^xt = m(4 + 3g— \/bq^ + lOg + 4). Finally, if the inner disk radius is inside the 
interval (2m, 2.34m), the naked singularity is of class I with q G (—1, —0.5527). In the case 
of class III and II naked singularities, it is possible to find the exact value of the quadrupole 
parameter by measuring the inner radius of the accretion disk. This is not true in the case 
of class I singularities, because all the circular orbits are stable in this region. However, if we 
could measure the radius and the angular momentum of a particular circular orbit located 
inside the region (2m, 2.34m), the quadrupole parameter could be determined by using the 
expression (129|) for the angular momentum. In this manner, we see that it is always possible 
to determine the quadrupole parameter from the physical properties of the accretion disk. 
The fact that in the gravitational field of naked singularities there exist circular orbits near 
the outer singularity is explained by assuming the presence of repulsive gravity. 

The study of radial geodesics starting at an arbitrary polar angle 6 shows that the 
quadrupole induces a deviation from the original radial direction in such a way that the 
test particles tends either towards the equator or to the axis of symmetry, depending on the 
value of the quadrupole. We also calculated the proper time and the coordinate time along 
radial geodesics on the equatorial plane and established the analogies and differences with 
respect to radial geodesics in the Schwarzschild spacetime. 

Our results show that the presence of a quadrupole can drastically affect the motion of 
test particles in a such a way that it is possible to determine the quadrupole by studying 
the properties of the test particle trajectories. If we were to compare our results with 
observational data from astrophysical compact objects, we would immediately notice that 
an important astrophysical parameter is missing in our analysis, namely, the rotation. We 
expect to perform such an analysis in a future work by using a stationary generalization of 
the g—metric as, for instance, the one derived in 321. 
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Appendix: Calculation of the ADM mass 


Consider the space-time (A4,g) foliated by a family (Etjten of spacelike hypersurfaces. 


Let 




(A.l) 


where is the boundary of a domain "iC [d'f' is assumed to be a timelike hypersurface). 
The ADM mass of the slice S* is dehned by 


Madm = 


1 


lim 




(A.2) 


1 f)7r j 

where / is a flat background metric on h is the induced metric of the hypersurface D 
stands for the connection associated with the metric /, S'* stands for the components of unit 
normal to ■s/s(fy denotes the surface element induced by the spacetime metric on y't 
with Sab being the induced metric, y°- = are some coordinates on and s = det Sab- 

Let us take for Et the hypersurface of constant coordinate t. Then, in the case of the 
g—metric ([T]) we have for the induced metric on the hypersurface Et 

n \ -q-1 / 2 ■ 2 n \ -'?(2+g) 

^ 2m \ /, m sm (J\ 


hij = diag 


1 + 


— 2mr 


.ryi-2f)-''(l + ^) 


-<?(2+g) 


2m\ 9 




(A.3) 


whereas for the flat metric we obtain 


fij = diag 


1 2 2 • 2 /I 

1, r , r sin 6 


(A.4) 


For we take the sphere r = const on the hypersurface E^. Then y°' = {0,ip), ^/s = 
r^sin6* and S* = 5*. Accordingly, Eg. (lA.2p becomes 


A^adm = 


lim 


IGtT t—)-cxd / , 

•''r=const 


[D^hrj — Dr{f^'‘hki)]r^ svaOdOdif^ 


where 


and 


Drif^hkl) = {hrr + T + T 


dr 


-2 ■ -2 
sm 




D^h^j Df^hj-j 

(jr r sm 6 


(A.5) 


(A.6) 


(A.7) 
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Here we used the fact that is a scalar held. In obtainging Eq. flA.7p . we used the 

formula 



— Ah — 2r^ h 

(A.8) 

DohrO 

^ Or^OOi ^ 

(A.9) 

D h 

^ip’ ^rip 

= —h —V^h 

Q^'^vip ^ ipr'^W ^ iptp'^rr’) 

(A.IO) 


with the non-zero components of the Christoffel symbols of the connection D with respect 
to the coordinates x* given by 


Tee =-A r;^^ = -rsinV, = = = - cos0sin0, = cot 0. (A.ll) 

Accordingly, from Eq. flA.6l) and Eq. flA.71) . after a simple but tedious calculation, one 
obtains that 


D^hr, - - q{2 + q)m^ (^1 

where we have used the hrst-order approximation 


2m\ ^ ^ 
r / 


3 2mq\ 

^ J 


sm^0, 


(A.12) 


1 + 


sin^ 0 \ 
— 2mr J 


-g(2+(j) 


1 - 


q{2 + q)rn? sin^ 6 
— 2mr 


(A. 13) 


Thus, by substituting Eq. flA.12p into Eq. flA.Sp and calculating the limit after integration we 
obtain Madm = m(l -|- g). 
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